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Abstract 

We extend the planar Pfaffian formalism for the evaluation of the Ising par- 
tition function to lattices of high topological genus g. The 3D Ising model on 
a cubic lattice, where g is proportional to the number of sites, is discussed 
in detail. The expansion of the partition function is given in terms of 2 2g 
Pfaffians classified by the oriented homology cycles of the lattice, i.e. by its 
spin-structures. Correct counting is guaranteed by a signature term which de- 
pends on the topological intersection of the oriented cycles through a simple 
bilinear formula. The role of a gauge symmetry arising in the above expansion 
is discussed. 

The same formalism can be applied to the counting problem of perfect 

matchings over general lattices and provides a determinant expansion of the 

permanent of 0-1 matrices. 

PACS Numbers : 05.20 - 64.60 - 87.10 
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I. INTRODUCTION 



The evaluation of the matching polynomial of a general graph with weighted edges is at 
the same time a root problem for discrete mathematics, statistical mechanics and mathemat- 
ical chemistry. Even in its simplest version, the so called dimer covering problem, in which 
the sites of a lattices have to be covered by non-overlapping arrangements of dimers, the 
evaluation of the perfect matching polynomial is a fundamental problem of lattice statistics 
]T|-^|]. For planar graphs, e.g. 2D regular lattices, the counting problem is easily reduced 
via Kasteleyn's theorem on lattices orientation to the evaluation of a finite number of Pfaf- 
fians mil]. Such a computation requires a number of operations which is polynomial in 
the number of vertices and is considered to be a tractable problem. For instance, the exact 
analytical solution of the regular 2D Ising model [[jj can be easily obtained by expressing the 
high temperature loop counting problem in terms of a dimer covering generating function 
over a properly decorated lattice . The periodic nature of Kasteleyn's orientation 

allows for the evaluation of the associated Pfaffian by diagonalization. Similarly, the Pfaf- 
fian method has been used in mathematical chemistry [Hj to derive the asymptotic number 
of dimer coverings for any regular surface lattice. Such a number is strictly related to the 
efficiency of adsorption processes of dimer molecules over surfaces, or to the degeneracy of 
double bond arrangements in planar organic lattices (the so called Kekule structures). 

In the case of non-homogeneous planar lattices, though the closed-form analytical solu- 
tion is in general impossible to obtain, the Ising and the Dimer problems remain tractable 
in algorithmic sense [[Tl|1 . 

The nature of the matching problem changes completely if one considers non-planar 
graphs or lattices ||. In discrete mathematics, it is known that the counting problem 
becomes ^P-complete |12] and no exact polynomial algorithm exists for the enumeration of 
coverings. 

In statistical mechanics and mathematical chemistry, the interest in non-planar lattices 
hinges on the fact that they are equivalent to higher dimensional lattices. The 3D cubic 
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lattice can be considered as a handlebody 2D lattice of topological genus g = 1 + N/4 
where N is the number of sites. A non vanishing ratio g/N for iV — > oo is related to an 
effective dimension D > 2 of the lattice, at least as far as its computational complexity is 
concerned. No exact solution exists for any non-planar lattice model, the simplest case being 
two coupled 2D Ising models. Similarly, no exact evaluation of dimer coverings over non- 
planar lattices is available. Of course, there exist several powerful probabilistic algorithms 
and approximate theories which provide quite accurate information, however the issue of 
understanding the onset of intractability is a basic open one. 

In this work we give an explicit formalism which generalizes Kasteleyn's method to 
arbitrary non-planar graphs. A first step in this direction was obtained in ref. [13[] in which 



the complete solution for the Ising model on a highly symmetric finite lattice of genus g = 3 
and N = 168 vertices was presented. Here we shall extend such formalism to any lattice and 
provide a general algorithmic procedure for the 3D cubic lattice. Scope of the paper is to 
link the combinatorial Pfaflian representation used for planar lattices with the topological 
features of non planar lattices. As a result we find an expansion for the 3D partition function 
in which the role of spin variables is played by a smaller set of binary topological excitations 
describing spin structures of the embedding surface of the lattice. 

Already in 1963, Kasteleyn |]|| noticed that the matching polynomial and the Ising 
partition function could be written as a weighted sum of 2 2g Pfaffians. In particular, since 
that time it is known that each Pfaflian can be associated to an element of the group 

(z 2 y x (z 2 y. 

In what follows we show that the Ising partition function can be written as Z = 
(2 cosh(/3J)) 3Ar Zo(A), where J is the spin-spin interaction energy, X = tanh(/3J) is the 
activity of a bond at inverse temperature f3 and Z (X) is the dimer covering generating 
function given as a series of Pfaffians with a topological signature. The final formula we 
shall prove is 

{e fe =0,l} k=l 



where the variables {e& = 0, 1} encode the orientation of the 2g elementary homology cycles, 
I[u)ki ^k') is the topological intersection matrix of the homology cycles Uk and <E> represents 
the orientation of the lattice. 

The paper is organized as follows. In Sec. II we outline some basic results concerning the 
combinatorial approaches to the 2D Ising model and we briefly remind the main steps of 
the so called Pfaffian method. In Sec. Ill we give a thorough description of the topology of 
the 3D cubic lattice, thereby fixing the notation. Sec. IV is devoted to the generalization of 
Kasteleyn's theorem and to the description of the gauge symmetry that such a generaliza- 
tion introduces in the problem. In Sec.V we analyze the set of cycles an co-cycles in terms 
of which the partition function will be expressed. The construction of a topological inter- 
section formula which gives the sign of Pfaffians in the expansion of the partition function 
is given in Sec. VI. The final constructive procedure is then presented in Sec. VII. Sec. VIII 
contains some preliminary results on the Pfaffian expansion whereas in Sec. IX we discuss 
the application of the formalism to the Dimer Covering and the Permanent problems. 

Throughout the paper a few numerical results will be given in order to provide some (very 
preliminary) physical insight. The analysis of the physical consequences of the formalism 
together with the discussion of the technicalities involved will be the subject of another 
paper fl26f . 



Independently, in refs. |L8[ some general results that partially overlap with ours are 
proposed. 



II. REVIEW OF COMBINATORIAL METHODS 

Despite the fact that the original Onsager solution to the 2D Ising model relied on 
the algebraic Transfer Matrix Method 0, the combinatorial solutions which have followed 
provide a more direct geometrical insight into 2D critical phenomena and field theories. 

While the transfer matrix method can be defined in any dimension, the combinatorial 
approaches strongly depend on the topology of the space where the lattice is immersed. Very 
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schematically, in 2D the sum over spin configurations is recast as a sum over closed curves 
(loops). Such curves are endowed with both an intrinsic topology and with the extrinsic 
one of R 2 . Since the Ising action depends only on the extrinsic geometry of loops, one has 
to avoid double counting and a proper cancellation mechanism, a topological term, has to 
be introduced in the sum. Such an approach has been developed by Kac and Ward WA 



and provides probably the most natural way of taking the continuum limit toward a field 
theoretical analysis Hl5|Jl6f . 

In 3D, the generalization of the above method encounters enormous difficulties due to 
the variety of intrinsic topologies of surfaces immersed in 3D lattices. Despite the deep work 
done in the attempt of recasting the critical 3D Ising problem as a string theory ||17|| , the 
problem remains unsolved under many aspects. 

Here we generalize the 2D (planar) Pfaffian or Dimer Covering approach to the Ising 
model, a purely combinatorial and basic tool of discrete mathematics that has many appli- 
cations in counting problems |J. 

In 2D, this approach relies on the equivalence between loop counting and dimer cover- 
ings (also referred to as perfect matchings) over a suitably decorated lattice. Once such a 
relationship is established the Pfaffian methods turns out to be simple both for the deriva- 
tion of exact solutions (in the cases of periodic lattices) and for the definition of polynomial 
algorithms on 2D heterogeneous models P,[H|. 

Let us briefly remind how the method works in the 2D case. The interaction energy of 
the Ising model on a planar square lattice A 2 d is given by 

Ni N 2 Ni N 2 

h = —j\ ^2 ^2 a j,k a j,k+i — J2 ^2 X! (2) 

j=l k=l j=l k=l 

where N\,N 2 are the number of sites in the two orthogonal directions, J\, J2 are the spin- 
spin interaction energies and cr^ = ±1- The partition function Z = J2{a=±i} ex P( — flH) can 
be written as 



Z = (cosh(^Ji) cosh((3J 2 )) NlN2 Yl 

{a=±l} 



Ni N 2 

3=1 k=i 



X 



Ni N 2 



II IK 1 + X 2CTj,kVj,k+l) 



(3) 



j=l fe=l 



where = tanh(/3Jj) are called activities of the bonds. Expanding the product and evalu- 
ating the sum over {a = ±1}, all the terms containing odd powers of cr's give no contribution 
whereas all even powers may be replaced by 1. It follows that the partition function acquires 
a clear interpretation as generating functions of closed loops with p horizontal and q vertical 
bonds with no overlapping sides. In fact denoting with N pq the number of such loops, we 
have 



In turn, the above expansion can be mapped onto the problem of evaluating the gen- 
erating function of dimer coverings (the so called weighted matching polynomial) over a 
new "counting" lattice A^d obtained by substituting each site of the original lattice with a 
cluster of six sites (two triangles with a joining bond) and by assigning activity 1 to the new 
decorating bonds while retaining the activity of the bonds inherited by the original lattice. 
The (eight) possible configurations of loop bonds at any Ising site are in 1 — 1 correspon- 
dence with perfect dimer configurations on the decorating cluster. Therefore the sum in (|J) 
coincides with the generating functions of perfect matching over the decorated lattice. 

Finally, in order to compute Z we orient the lattice according to the Kasteleyn prescrip- 
tion by assigning arrows to each bond in such a way that for any closed circuit £ on Af D , the 
number of bonds of i oriented clockwise is of opposite parity to the number of sites enclosed 
by i. The Kasteleyn rules define completely the orientation for planar lattices, whereas for 
non-planar lattices , i.e. lattices which can be immersed on a surfaces of non-trivial topo- 
logical genus, we need further sign fixing for loops not homologically trivial (i.e. without an 
interior). The dimer covering generating function can then be expressed as a weighted sum 
of Pfaffians of the antisymmetric adjacency matrix with elements given by the activities of 
the bonds and signs determined by their orientation. In virtue of Cayley theorem, Pfaffians 
are computed as square roots of the determinant of such matrices. Thus, the Ising partition 



Z = (2 cosh(/5Ji) cosh(pj 2 )) NlN * J2 N P g X i x ! 



(4) 
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function can be written explicitly as a determinant which for uniform interaction energies 
can be further block diagonalized by Fourier Transform. The final calculation of a 6 by 6 
determinant leads to the exact closed form expression of the 2D Ising partition function. A 
thorough discussion of the above procedure can be found in ||. 

Below we shall concentrate on the generalization of the above construction to the cubic 
3D lattice. The procedure is however general and can be straightforwardly generalized to 
any non planar lattice. A first explicit example was presented in for the case of group 



lattices with non-trivial topological genus. The same inductive reasoning used in ref. f!3 
leads to a simple topological expression for the coefficients in the Pfaffian expansion. 



III. THE 3D CUBIC LATTICE AND EMBEDDING SURFACE 

We consider 3D cubic lattices A of sides A 1; A 2 , N 3 with N = NiN 2 N 3 sites and peri- 
odic boundary conditions. Each vertex V is identified by a triple of periodic coordinates 
{ni,n 2 ,n 3 }, m = O...Aj — 1 with V(rii,n 2 ,n 3 ) = V(mod(ni, Ni), mod(n 2 , A 2 ), mod(n 3 , A 3 )). 
The sites can also be labeled in sequential order by the single index q = q(ni, n 2 , n 3 ) = 
mod(n 1 ,A 1 ) + A!Hiod(n 2 ,A 2 ) + A!A 2 mod(n 3 , A 3 ) with the inverse relations, n\ = 
mod^AO, n 2 = mod(^,A 2 ) and n 3 = mod( g ~ n ^ 2jVl , A 3 ). In what follows nota- 
tions and operations over the integers %, n 2 ,n 3 have to be understood modulo Ni,N 2 ,N 3 
respectively. 

The lattice A is invariant under translations Dj : rij — > rij + 1. 

The set of N b = 3N bonds Li(q), i = 1,2,3 of A connects couples of neighboring sites 
{V(q), V(Diq)}, thus defining the adjacency or incidence matrix A of A, A q<q ' — 1 if q and 
q' are connected by a bond and A q A > = otherwise. 

We call plaquette a square face (ni, n 2 , n 3 ) = F^q) of A identified by the sequence of 
vertices, V(q), V(D i2 q), V(D i2 D i3 q), V(D i3 q), where with the notation ii,i 2 ,i 3 we denote a 
generic cyclic permutation of the indices 1,2,3. A contains 3 classes of N plaquettes Fi(q), 
orthogonal the axes i, {i = 1, 2, 3). 
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The parity of a site is given by p(q) = (— i) n i+"2+"3. A is a bipartite lattice in that edges 
always connect vertices of opposite parity. 

In order to implement the dimer method we construct an orientable surface E without 
boundary, which contains all the sites and bonds of A and is the union of a subset of square 
plaquettes of A. The number Nf of such plaquettes is Nf = N b /2 = 3N/2, each bond 
belonging to two plaquettes of the surface and each plaquette containing four bonds. It 
follows that A" and at least one of the numbers Aq, N2, A3 need to be even. For simplicity 
we shall assume N { = 2M if so that N = 8M with M = MiM 2 M 3 . 

As we shall see, all the above conditions can be matched by a definition of E which 
preserves part of the symmetries of the original lattice. 

The topological genus g of the surface, evaluated by Euler's formula, is N — N b + Nf = 
2(1 - g), from which it follows 5- = l + 2M = l + N/4. 

The definition of E requires that a plaquette (q) belongs to E only if rii 2 + rii z is odd, 
and we shall call faces such plaquettes. The final result for E consists in a square beams 
periodic structure shown in Fig. 1. 



In order to proceed in the generalization of the Pfaffian method it is useful to recall some 



Sites, bonds and faces of E generate Abelian groups, additive modulo 2, of non-oriented 
chains Cfc(E,Z 2 ) of dimension k = 0,1,2 respectively. For any c G C^(E, Z 2 ), we have 
c + c = , where is the identity of the group. 

The linear boundary operator 5 maps chains of different dimensions 5 : C&(E,Z 2 ) — > 
Cfe_i(E, Z 2 ) and is defined as follows 



A. Combinatorial Topology of the 3D cubic lattice 



basic notions of combinatorial topology p0| . 



SV(q) 



1 



SLi(q) 



V(q) + V(D iq ) 



SF h (q) 



L i2 (q) + L i3 (D i2 q) + L i% {D i& q) + L h {q) 



(5) 



S 



Clearly, 5 2 = 0. A chain Ck in the kernel of 5, i.e. such that 5ck = 0, is said to be closed 
. A circuit is a closed sequence of vertices connected by edges, each edge is an element of 
Ci(E, Z 2 ), the sum of edges of a circuit is a closed 1-chain c\. 

A chain Ck-i = Sck G Cfc(E,Z 2 ) is called boundary. Since 5 2 = 0, a boundary ci is 
necessarily closed and reduces to a sum of circuits. 

We will keep the definition of SFi(q) as given by (fjD also in the improper case in which 
Fi(q) is not a face and hence 8Fi(q) is not a boundary but only a circuit on E. 

Boundaries and closed chains of Cfc(E,Z 2 ) generate the subgroups -Bfc(E,Z 2 ) C 
Z fc (E,Z 2 ) C C fc (E, Z 2 ), respectively. 

The homology group is thus defined by H k (T,, Z 2 ) = Z fc (E, Z 2 )/_B fe (E, Z 2 ) so that there 
exists a projection n from closed chains to the elements of Hi, 7r : Zi(E, Z 2 ) — > Hi(E, Z 2 ). 
Boundaries are topologically trivial since they are mapped by it onto the identity of the 
homology groups. The elements of i?i(E, Z 2 ) are called cycles and the homology group is 
generated by a base of 2g elementary cycles which are equivalence classes of closed chains 
under the addition of boundaries. A closed chain Ck G Z fc (E, Z 2 ) such that nek — jk G 
i/jfc(E, Z 2 ) is called a representative of 7^. 

The multiplicative functionals on the elements of Cfc(E,Z 2 ) and iffc(E,Z 2 ) with values 
±1 constitute groups of co-chains and co-cycles denoted by C fc (E, Z 2 ) and if fc (E, Z 2 ) respec- 
tively. The definitions of i/fc(S, Z 2 ) and if fe (E, Z 2 ) are independent from the triangulation 
of the surface E and in particular will be valid for the decorated lattice T defined in the 
following section. The properties of the groups H k (T l ,Z 2 ) and H k (T>,Z 2 ) depend only on 
the topological features of E and not on the choice of the tessellating lattice. 

The symmetry properties of E can be seen more clearly by considering its embedding 
into the subset or box 5 of R 3 defined by0<a;<iV 1 ,0<?/<iV2,0<,2<iV3 with 
periodic boundary conditions. The complement H — E is the union of two open and disjoint 
subsets S + , H_ congruent under the translation S = D1D2D3 and invariant under = Df. 
S and {Tj} generate the symmetry group of E while Di alone is not a symmetry of E. The 
closure of E+, H_ is given by S + = S + U E, S~ = S_ U E and H + R S~ = E. We may 
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conventionally regard S + as the interior of E. On S + , the bond Li(q) is convex/ concave 
depending on whether rii(q) is odd/even, hence, moving along S + with Di, we encounter 
alternatively convex and concave bonds. On H + , each face has 2 concave and 2 convex bonds 
which interchange if we focus on 5~ rather than S + . 

IV. KASTELEYN'S ORIENTATION AND GAUGE SYMMETRY 
A. Decorated counting lattice 

In order to write the Ising partition function as a dimer covering generating function we 
first construct a decoration T of E obtained, following Fisher's prescription ||, by replacing 
each site of coordination q with a graph of 3(q — 2) points and q — 2 triangular faces, as 
shown in Fig. 2 (for the case of interest here q — 6). In our case, the decorated surface F 
contains 12iV vertices. 

Each face of E maps into a face of T composed of 4 bonds inherited from E plus 8 bonds 
introduced by the decoration. Each site of E originates 4 triangular faces in T, Fig. 2B. 
As we shall see, the above decoration is also equivalent (with respect to dimer covering 
configurations) to a locally non planar one in which each site of E yields a complete graph 
of 6 vertices, Fig. 2C. 

Clearly, the homology groups of T and E coincide and will be identified. 

B. Proper orientations 

The orientation over bonds is an additional geometrical structure $ defined by a function 
<pj(q) = ±1, j = —3,-2,-1,1,2,3 such that <fij(q) = 1 (resp. —1) corresponds to a bond 
oriented from V(q) to V(Djq) (resp. from V(Djq) to V(q)), and (f)-j(q) = —Dj 1 (pj(q). 

An orientation of the bonds of T or E is said to be proper if it satisfies the condition of 
Kasteleyn's theorem for planar lattices, i.e. if by moving anticlockwise along the perimeter of 
a face we encounter an odd number of oppositely oriented bonds. Unless otherwise specified 



10 



we assume in the following that all orientations $ are proper. Consider now a boundary on 
T or S consisting of a single circuit and containing in the interior k points. By induction 
it follows that moving along it anticlockwise we encounter a number of oppositely oriented 
bonds with parity opposite to that of k. A proper orientation of E defines a proper one for 
T as follows. Bonds of T inherited from E are given the same (f)j(q). Bonds forming the 
boundary of an ornating triangle are then all clockwise oriented so they appear anticlockwise 
in the adjacent faces of T inherited from E and lead to a proper orientation. Sites of E and 
corresponding decorating clusters of T fall into 2 subsets of even and odd sites according to 
their parity p(q). 



C. Gauge group 

Given a proper orientation $ of T, it is possible to derive different but equivalent ones by 
reversing all the bonds of T which are incident to a given site V(q), according to the gauge 
operation (see Fig. 3) 

g(q)--<f>j(q) -> -<M<?) , Y? (6) 

The operations {g(q)} generate the gauge group Q. The orientations generated by Q are 
equivalent and will be identified. By using Q we can fix arrows on T such that they relate 
even and odd clusters by a mirror reflection and reversal of orientations of all the bonds. 



D. Generalized Kasteleyn rule for dimer coverings 

A dimer covering of T is defined by a 2-coloring (say black or white) assignment to the 
edges of T such that each site of T belongs exactly to one black bond. By superposing two 
dimer coverings we obtain a set of closed black circuits which cover all vertices of T. By 
moving around a circuit we encounter an alternating sequence of black bonds from each of 
the 2 coverings, therefore the length of the circuit must be even. If the circuit is a boundary 
it must enclose points which are themselves connected by bonds forming black circuits, thus 
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the number of points inside must also be even and by moving around the circuit in whatever 
direction we encounter an odd number of opposite arrows. 

Dimer configurations are in one to one correspondence with terms of the Pfaflian of an 
antisymmetric incidence matrix M($,X) which will be uniquely defined further on by the 
pair (r, $). In what follows we shall use 1-chains composed by black circuits only, origi- 
nating from the superposition of dimer coverings. For simplicity, since no confusion may 
arise, we use for the black subgroups and the corresponding dual groups the same nota- 
tion .Bi(E, Z 2 ), Zi(E, Z 2 ) and i?i(E, Z 2 ) = -Zi(S, Z 2 )/5i(S, Z 2 ) introduced in the previous 
section for the full groups. 

Adding a boundary to a black circuit c does not alter its orientation. The latter depends 
therefore only on the image 7 = tt(c) G -Hi(E, Z 2 ) and defines a functional $(c). From the 
bond orientation function <pi(q) we define the orientation function $j 1 (q r ) for the boundary 
of the faces SF^ (q) as the multiplicative and gauge invariant functional 

= -(t> i2 (q)(l>i 3 (D i2 q)(/)_i 2 (D i2 D i3 q)(t)_ i3 (D i3 q) . (7) 

For a black boundary, Kasteleyn's rule implies 

= 1, i = l..3,q = l...N . (8) 

Given an anticlockwise sequence of bonds L^.^qj) forming the circuit c\ G Ci(E, Z 2 ) and 
connecting the points V(qj), j = 1, ...m, the orientation $(cj.) over i?i(E, Z 2 ) is defined by 

$( Cl ) = -n- =1 fe .fe) . (9) 

On a lattice of non trivial genus g (g > 0) such a functional is not an element of 
if 1 (E, Z 2 ), i.e. in general $(ci + c 2 ) 7^ $(ci)$(c 2 ) , as it can be readily verified in the simple 
2D Ising lattice with periodic boundary conditions but also on other finite lattices ||13| . The 
argument runs as follows. 

The single black circuit c li2 given by the sum c\ + c 2 of two black circuits c\ and c 2 of 
parity $(ci) and $(c 2 ) respectively, which intersect p times, has parity 
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$(c 1)2 ) = (-l) p $( Cl )$(c 2 ) . (10) 

This fact will turn out to be essential in determining the signs with which the Pfaffians 
have to appear in the expansion of the generating function. Let us consider the specific 
example of the sum of two black circuits C\ and c 2 of parity $(ci) and $(c 2 ) with represen- 
tatives a + i for Ci, b + 1 for c 2 where i is a common bond and a, b are paths closed by £. A 
representative for the product c i2 is now a + b. 

We move along c\ in a specific direction encountering n a opposite arrows along a and 
r = 0, 1 opposite arrow coming from I with a total of n a + r. Similarly we proceed along c 2 
in such a way that when we move along ci i2 , a and b are run in the same direction and meet 
a number + 1 — r of opposite arrows since i is run in the opposite direction in c 2 with 
respect to c\. The total number of opposite arrows along c 1)2 is now n a + n b while the total 
on ci, c 2 is n a + n& + 1. This result is equivalent to (|10|) for p = 1 and can be generalized to 
multiple intersections. 



E. Gauge tree, spin structures and Pfaffian matrices 

The gauge symmetry G can be used to fix the orientation on a subset of bonds forming a 
spanning tree T in T, see Fig. 4. It will be seen that both the gauge choice and Kasteleyn's 
rules lead to a number of orientations which is 2 2g . Such orientations break all symmetries 
of S and are in one to one correspondence with the 2 29 spin structures of a surface of genus 
9- 

For each element of the base of iYi(£,Z2) there exist two possible orientations, for a 
total of 2 2g different global orientations of S. Given $ for the cycles of the homology base, 



it is possible to derive the parity of any circuit by the intersection formalism related to flTop. 

$ defines (modulo gauge transformations), an antisymmetric matrix M($,X) of dimen- 
sion 12N x 12N with elements labeled by the sites of T. Changing the gauge produces a 
new matrix M'($,X) = KM($>, X)K~ 1 , with K diagonal, which has the same Pfaffian of 
M($). 
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If q,q' are not connected we set M M /($,X) = 0. If q,q' are connected by a decorating 
bond Li(q) , we have M M /($, X) = 4>i(q) = ±1 depending on the orientation of the bond. If 
q,q' are connected by a bond inherited from A, we set M gi9 /($, X) = <fri(q)X = ±X where 
X is the so called activity of the bond, X = tanh(/3J) with (3 inverse temperature and J 
spin interaction energy. The diagonal 12 x 12 blocks of M($, X) describe the decorating 
clusters displayed in Fig. 2B and are given by 
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(11) 



We have thus defined all the elements of M($,X) which can be written as 

M($,X) = A + X 



(12) 



where A is block diagonal consisting of N/2 fixed blocks as in ((TTJ) corresponding to even 
sites and N/2 opposite blocks corresponding to odd sites. -B($) contains the <f>i(q). For 
computational purposes the matrix M($,X) can be replaced by an equivalent one of di- 
mensions 12M x 12M by a folding procedure which we can sketch as follows. In this section 
I n denotes the identity matrix of rank n. 

i) The 12 AT sites are grouped into 2 subsets each of 6N sites . The first subset Ext 
includes the internal sites of each of the N decorating clusters, which label in ( |Tl~|) rows 
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and columns 6, 12. Each of these rows and columns have 3 non zero entries. The second 
subset Int contains sites with 2 decorating bonds (the 2 non zero entries in ( PD ) and one 
inherited from E which appears in B. Schematically M ($, X) can be written in blocks 
of the form 



M($,X) 



Mu Mi 2 
M 21 M 22 



(13) 



where i,k = 1,2 label the subsets Int respectively. But now Mu is a block diagonal 
matrix with Det(Mn) = 1 and which can be easily inverted. We write then 



det(M) = det(Mn) det(M 22 - M 21 Mf 1 1 M 12 ) = det(M a ) 



(14) 



M a = M 22 - M 21 Mil 1 M 12 is now a QN x QN matrix where the decoration now reduces 
to complete graphs of order 6 shown in Fig. 2C. 

ii) The 6N Ext entries in M a are partitioned again into 2 subsets labeling even and 
odd sites respectively thus exploiting the fact that £ is a bipartite graph. Schematically M a 
can be written as 



M, 



M EE M EO 
Mqe Mqo 



(15) 



Now Mee-, Moo do not contain X since an inherited bond always connects sites of opposite 
parity while Meo-, Moe are instead linear in X. We may apply again (|I~4"D and obtain a 
3N x 3iV matrix of the form 



M b ($,X 2 ) = A b + X 2 B b ($) 



(16) 



such that det(M a ) = det(M 6 ). 

It can be verified that the diagonal matrix Q — Ijv ® Diag(l, 1,1,— 1,-1,— 1) obeys the 
relations: 



QA c Vl = A^ 1 , VlB c n = B^ 1 , 



(17) 
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and that A b ,B b share the same spectrum with eigenvalues ±z, ±i(2 — y3) > ±i(2 + \/3) 
independently of $ and det = det Q = 1. The essential information is now contained in 
B b and X since Aj, is a fixed numerical matrix. 

By Cayley's theorem the determinant of M($, X) or M&(<&, X 2 ) is the square of an even 
polynomial Pf($,X) in the activity X called Pfaffian of M($,X). 

As we shall see, the ambiguity of sign in the extraction of the square root can be solved 
by imposing all Pfaffians to be = 1 in the high temperature limit X —>■ 0. However the 
polynomial Pf($,X) may have zeros on the real axis in the physical range X = 0, 1 and 



may change sign when reaching X = 1. By generalizing the results of ref. |13| we find that 
the sign in X = 1 is given by the function <r(a) defined below in (^3|), having a precise 
topological meaning, 
iii) We have : 

det M 6 ($, X 2 ) = det(A b + X 2 B b ) = det(l 3N + X 2 A b x B b ) = det(l 3N + X 2 U b ) (18) 

where U b = A b ~ 1 B b . 

Let 6 = QA b so that 2 = l3jv and 

QU b & =®A b : 1 B b Q = VLB b VLA b = B b - 1 A b = U b l (19) 

is again a block 6x6 matrix which can be explicitly diagonalized so that by changing 
basis both B and U can be written in the following block form where briefly I = Ium — I3A72 

(20) 





I 




Q R 


= 




,u b = 




-I 




S T 



From (|19"1) we see that Q,R, S,T satisfy the identities 

Q 2 — RS = I, Q = RTR-\ Q = S~ l TS . 

We have now 



(21) 



det (I + X 2 U b 



det 



I + X 2 Q X 2 R 
X 2 S I + X 2 T 



det 



I 

S 



I+X 2 Q X 2 RS 

X 2 l+X 2 S- l TS 



I 

s- 1 



) 



(22) 
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but now S can be included in the changes of basis and we can assume S = 1, T = Q, R = 
Q 2 — I. Therefore we may write simply 

-(X 4 + 1)I-2X 2 Q 

1 I+X 2 Q 

= det ((X 4 + 1)1 + 2X 2 Q) = X 24M det ((X 2 + X~ 2 )I + 2Q) (23) 

where we subtracted from the upper block row the lower block row multiplied by Q + X~ 2 . 
In the last form the determinant is evaluated for a matrix (X 4 + 1)1 + 2X 2 Q of rank 12M, a 
factor 8 down from M($,X), with a considerable gain in the computational speed. Clearly 
Pf($,X)X~ 12M is invariant under the map X — > X -1 . 

V. CYCLES AND CO-CYCLES OVER £ 

The black elements of Z 2 ) are generated by two classes of elementary cycles 

E il (q),Oi 1 (q), referred to as even or odd cycles forming the sets E,0. Such cycles are 
homologically equivalent over T to cycles of the form 

L i2 + L i3 (D i2 q) + L i2 (D i3 q) + L i3 (q) , (24) 

where n^n^ are even/odd for cycles G Ej O 

Cycles G {E,0} are thus generated by the boundaries 5Fi(q) of the 12M = 3X/2 
plaquettes of F which are not faces of S (see Fig.l). 

However, these cycles are not independent since their number is 6M > 2(1 + 2M) = 2g. 
In order to perform a correct counting we consider the following product of 6 plaquettes for 
n\ + 1, n 2 , n 3 even (remind q = (ni, n 2 , n 3 )) 

E 1 (D 1 q) + E^q) + 5F 2 (q) + 5F 2 (D 2 q) + 5F 3 (q) + 8F 3 (D 3 q) , (25) 

i.e. the boundaries of the plaquettes of the elementary cube C(q) obtained by translating 
of the site V(q) by one unit along the positive direction of all three axes. By expand- 
ing the product (p5|) each bond appears twice, whence follows the equivalence modulo a 
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det(I + X 2 U b 



det 



I+X 2 Q X 2 (Q 2 - 1) 
X 2 I+X 2 Q 



det 



boundary of Ei(Diq) and Ei(q), which are identified in E. Indeed the parity of the last 
four plaquettes indicates that they must be faces of E and are topologically trivial. Both in 
F 3 (ni,n 2 ,n 3 ),F 3 (ni,n 2 ,n 3 + l) the sum ni+n 2 is odd while in F 2 (ni, n 2 , n 3 ), F 2 (ni, n 2 + l, n 3 ) 
the sum ni + n 3 is odd. It follows that they have boundaries equivalent to the identity of 
H 1 (H,Z 2 ) and that 

Ei(Diq) = Ei(q) , m odd 

Oi(Diq) = Oi(q) , rii even . (26) 

Each cube of A gives an identity among cycles, but their actual role depends on the 
parity of ni, n 2 , n 3 . 

If ni,n 2 ,n 3 have all the same parity, the boundary of the cube does not contain faces 
and leads to the identity 

j2(E i (q) + E i (D i q)) = , (27) 

which connects 6 cycles of the same parity and must be added to the previous identities. We 
shall denote with £ and O the set of values of q restricted to ni, n 2 , n 3 all even and all odd 
respectively. If q e £, the cubes C(q) and C(Tjg) are contiguous and contain the common 
cycle Ei(Diq) = E^Tiq). It follows that the M identities deriving from the even cubes are 
not independent because in their sum each cycle appears twice. The same reasoning applies 
to the odd cubes. Thus we may write 

5>i(9) = 0, ]TO,(g) = . (28) 

Another class of identities connects even with odd cycles. Consider the plane of A with 
n 3 constant. Such a plane is tessellated by plaquettes which are faces if ri\ + n 2 is odd and 
otherwise have boundaries (on A but not on S) which are cycles G E 3 (q) or 3 (g). In virtue 
of the periodicity of A the plane has toroidal topology and the sum of the boundaries of all 
the plaquettes is 0. We thus have the identities 

e(n 3 ) = E 3 (n 1 ,n 2 ,n 3 ) = o(n 3 ) = J2 3 (n u n 2 ,n 3 ) . (29) 

ni,ri2even ni,n2odd 
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Replacing in the above expressions n 3 with n 3 — 1 we obtain the same identity due to 

Such a result can be further extended as follows. If n 3 is even, we have e(n 3 ) = e(n 3 — 1). 
Let's consider the product of the boundaries of the cubes C(ni, n 2 , n 3 ) where n 3 is a fixed 
even index and rii,n 2 are even 

[Ei(ni,n 2 ,n 3 ) + Ei(Di(ni,n2,n 3 ))} . (30) 

m,ri2 even, i=l, 2, 3 

Factors of the type E\{n\, n 2 , n 3 ), E 2 (rii,n 2 ,n 3 ) always appear twice in contiguous cubes 
and the products can be written as 

53 [E 3 (n 1 ,n 2 ,n 3 )+E 3 (n 1 ,n 2 ,n 3 + l)}=e(n 3 )e(n 3 + l) , (31) 

ni,n2even 

from which we have e(n 3 ) = e(n 3 + 1). It follows that e(n 3 ) does not depend on n 3 . The 
same result holds for all directions and for cycles of odd arguments. Form ( p9| ) we have 
e(rij) = o{rii) = I i: providing three exceptional identities Jj, i = 1..3 which connect cycles of 
different parity. 

Finally we have three exceptional cycles Ri(q) i = 1..3 given by 

R t (q) = L i( D il) ■ (32) 

n = 0,l,...,JVj-l 

We set Ri = Ri(0). Each Ri(q) can be expressed in terms of Ri and of E,0 cycles. Note 
that Ri(q) does not depend on rii and hence can be written in terms of two site labels and 
a direction, e.g. R\(q) = Ri(n 2 ,n 3 ). Geometrically Ri(q) is a straight circuit which winds 
around A by exploiting the periodicity, its length is N\ = 2M\ and is even (see Fig. 1). 
There is no ambiguity in extending formula (Q) to black Ri(q). The existence of Ri(q) follows 
from the periodicity conditions imposed on A and also exists in the 2D Ising model where 
they produce the analog of Ri(q) for i = 1,2. There is however no analog in 2D of the E, O 
cycles. 

We are now in position to check that the overall number of independent cycles is in fact 
2g. To each one of the M even/odd cubes corresponds one identity but only M — 1 of them 
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are independent. The 3M even/odd cycles we started with reduce to 3M — (M — 1) = 
2M + 1 = g independent ones, giving exactly an overall number of 2g. The three cycles 
that are eliminated by the identities are replaced by the exceptional cycles Ri so that 
the overall counting remains unaltered. (In the thermodynamic limit N — > oo, we expect a 
negligible contribution from the surface terms Ri). 

We can restate the counting problem by noticing that the total number of defined cycles 
amounts to 3M + 3M + 3 = 6M + 3 = 3g. However the total number of identities is given 
by 2(M — 1) + 3 = g and hence the number of independent cycles is once more 4M + 2 = 2g 
as expected. 

The generalized dimer method amounts to writing the partition function as a sum of 2 2g 
G-invariant Pfaflians associated to all possible different orientations $. 

We call fundamental orientation <&p that for which $(a) = 1 for all a £.ffi(£, Z2). 
This can be obtained by setting 

(j>i(ni,n 2 ,n 3 ) = (-1)** . (33) 

While the fundamental orientation makes no distinction among plaquettes, the antifun- 
damental §a requires <&i(q) = —1 for all E,0 cycles but $(i?i) = 1. There is no simple 
recipe for $^ analogous to (|3~3D since the way it appears depends on the spanning tree T. 

VI. TOPOLOGICAL INTERSECTIONS AND SIGN FUNCTIONALS 

A. Intersection of cycles 

The bilinear symmetric functional / [a, b] over the cycles a, b is given by Mod(p, 2) where 
p is the number of intersections of a, b over E. In our case a detailed analysis shows that 
cycles G E, O of same parity or same normal do not intersect. The general formula can be 
deduced from 

I[E 1 (m 1 , m 2 , m 3 ), 2 (^i, n 2 , n 3 )] = 5 ni , mi _i ^m2,n2- 1^7i3,m3+l j 

(34) 
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by imposing invariance under cyclic permutation of the axes and under the S, 7* symmetry 
operations. Two adjacent cycles do not intersect. 

The cycles Ri do not obey this rule and the intersections with the exceptional cycles are 
given by 

I[R k ,E l (q)} = iii^k 
I[R k ,O l (q)]=0 iii^k 

I[Ri, Ei(ni,ri2, n 3 )] = 5i,i(5n 2 ,o^n 3 ,o + ^2,^2-1^3,^3-1) 
I[Ri, Ei(ni, n 2 , n 3 )] = ^,2(^1,0^3,0 + ^1,^1-1^3,^3-1) 

I[Ri, E i (n 1 ,n 2 , n 3 )\ = 5 ij3 (5 nu0 5 n2 , + ^1,^1-1^2,^2-1) , (35) 

together with the expressions obtained by interchanging E with O. 

The definition of J [a, b] can be extended by linearity mod 2 to arbitrary pairs of cycles 
a, b. 

B. Topological excitation and signature 

The intersection formalism leads naturally to the notion of elementary topological exci- 
tations {ri(q)}, a set of minimal operators which change locally a proper orientation into a 
new one, still proper but inequivalent. r acts on the inherited bonds only and we can define 
r directly over E. Any two inequivalent proper orientations are connected by a sequence of 
r operations and hence {ri(q)} generate all orientations of E. In particular we may reach 
any orientation by repeatedly applying 7$ (5) to $ - 111 the following we use the equivalent 
notations $j(ni,n 2 ,n 3 ) = &i(q) = $(a) where a = Ei(q) or Oi(q). The elementary topolog- 
ical excitations are in one to one correspondence with the E, O cycles and we may adopt the 
same set of indices and distinguish between even or odd excitations, Ti(q) — > Tf(q),Tf > (q). 

The action of Ti(ni,n2,ns) = T^{n\ — 1,712,713) is given by 

0i(rai-l,n 2 + e,n 3 + f) -> -0i(ni-l,ra2 + e,n 3 + f) e,f = 0,l. (36) 
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The definition extends in an obvious way to the other directions and, in virtue of the 
invariance under the S, T{ operations, ( p6|) hold also for odd operators. Clearly, the values 
of (f> on bonds not appearing in ( [36]) are unaffected by rf(q). 

Given an elementary cycle, e.g. Ei(q), there exist four other intersecting cycles, e.g. 

2 {D^ l D 2 D^q) , 2 (Di 1 D 2 D 3 q) , 

Os(D^D^D 3 q) , 3 (D^D 2 D 3 q) . (37) 

Consider now the orientation $2(^1 — 1, n 2 + 1, n z ~ 1) corresponding to 02(721 — l,n 2 + 
l,n 3 -l) = 2 (D^D 2 D 3 1 q), 

$2(^1 - l,n 2 + l,n 3 - 1) = «^i(ni - l,n 2 + l,n 3 - l)0 3 (ni,n 2 + l,n 3 - 1) x 

0-i(ni,n 2 + l,7i 3 )0_ 3 (ni - l,ra 2 + l,n 3 ) • (38) 

The last factor changes sign under the action of T^(q) and hence the orientation of 
$2(^1 — 1, ^2 + 1, n 3 — 1) also changes and the same result holds for the other cycles of the 
example (|3"7|). 



The action of r i ' (q) changes $ into an inequivalent orientation which differs from $ 
only in the orientations of the local cycles intersecting Ei(q),Oi(q). The orientation of all 
faces of £ remain unchanged under r { ' (q) so that the Kasteleyn's conditions are always 
fulfilled. Each cycle can be given an active role if identified with 7$ ' (q) or a passive 
one if considered as a cycle changing parity under the action of r i ' (q). The functional 
J a G H 1 ^, Z 2 ) : #!(£, Z 2 ) -> Z 2 is defined by 

J a (b) = , a, b G Fx(S,Z 2 ). (39) 



C. Axial gauge 

We fix the gauge for the orientations of bond by selecting a subset of iV — 1 bonds forming 
a spanning tree T of A, that contains all the sites of A. A convenient gauge fixing is the 
Axial Gauge (as displayed in Fig. 4) 
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3 (O,O,n 3 ) = l n 3 = 0,1, ...,#3 -2 

2 (O,n 2 ,n 3 ) = (-1) 713 , /la = 0,1, ...,#2 -2, n 3 = 0, 1, iV 3 - 1 

0i(7ii,n 2 ,n 3 ) = (-l) ni , 7ii = 0,l,...,iVi-2, n 2 = 0, 1, N 2 - 1, n 3 = 0, 1, iV 3 . (40) 

The axial gauge leaves undetermined the orientations of 3iV— (N— 1) = 2iV+l = 16M+1 
bonds not belonging to T . However, the 12M — 1 Kasteleyn's conditions reduce the number 
of such independent orientations to 16M + 1 — 12M + 1 = 4M + 2 = 2g, as expected. 

From the recursive relation equivalent to (|10|) 

$(a + b) = $(a)$(b)(-l) /[a ' b] , (41) 

we deduce then the value of $ on all the elements of i?i(S, Z 2 ) starting from <&i(q). From 
(PH|) we see that <&i(q) is invariant under Q and the same is therefore true for all $(a). The 
set of all the <&i(q) determines therefore the global orientation $ (mod Q) of E. 

VII. GENERAL PROCEDURE AND PFAFFIAN EXPANSION OVER E 

The geometrical structure on which we define Pfaffians admits an alternative equivalent 
definition. We consider 2 cubic lattices A^, Ao each having M vertices in 1—1 correspondence 
with even/odd cubes C(n\, n.2, ti 3 ) where n\ + n 2 + n 3 = even/odd, respectively. The bonds 
of the lattices correspond to the even/odd E, O cycles. 

In place of spins we have orientation parity of the E, O cycles, satisfying the identities 
(p6|-p9D. Topology plays a marginal role in the classical 2D Ising lattice where the genus 
g = 1 leads to a sum over 4 pfaffians only. 

Here the reduction of the original sum on N = A(g — 1) spins to one on 2g functionals 
does not solve the problem but reduces the complexity of the task. The orientations on the 
sublattices A E ,A Q are not independent. 

We sketch now the general algorithm which fixes all gauges and gives all bond orientations 
4>i(q) in terms of a subset of 2g independent <p's by means of the Kasteleyn's conditions. 
We consider anly the homogeneous case where the activities of bonds are the same in all 

23 



direction, however the results can be straightforwardly generalized to any non homogeneous 
distribution of bond interaction energies. 

The actual steps are quite complex and can be summarized as follows: 
(1.) Compile a list of all 3iV <pi(q) considered as independent binary variables. 
(2.) Impose the axial gauge ([5]) and obtain a sublist of 2N + 1 terms 4>i(q) only. 
(3.) Impose the Kasteleyn conditions over faces. This can be achieved by solving iteratively 
(§) thus reducing the number of independent <pi{q) to 2g = 2 + N/2 forming a basis B for 
the 2 2a functionals in H 1 (E,Z2). From B we define by group multiplication the generic 
functional $(a) G Z 2 ), Va e#i(£, Z 2 ). 

(4.) Select a complete basis Q = {uo k ,k = 1..2g} of independent cycles out of 
E,0, Ri, R 2 , R3. A generic cycle can then be written as 

29 

a = E , e k = 0, 1 (42) 

k=l 

(5.) Define the self intersection function i2i(£, Z 2 ) -> Z 2 , see (0), ([35]), (|TD 

tr(a) = (-l)ELE^Vk^]e t e t , ( 43 ) 

(6.) The function X a in (^) defines an invertible duality map Z 2 ) — > H 1 ^, Z 2 ). The 

map cro X" 1 lifts the self intersection to Z 2 ). In any case we do not need to compute 

explicitly the inverse of X a . 



(7.) Given a as in (|42|), we compute X a and expand it in terms of the restricted basis B 
defined in step 3 , thus determining the orientations $(a), <pi{q) and all matrix elements of 
M($(a)) explicitly as functions of ex, ..■,e 2g . 

(8.) The dimer generating function Zq(X) is then given by the sum over Pfaffians 

Zo(I)4 E <7(a)Pf($(a),X) 

4 E (-^E^EXVK^le^, Pf ($(^ efc ^),X). (44) 

{e fc =0,l} fc=l 

As in eq.(|), the Ising partition function is simply given by 
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Z= (2cosh(/?J)) 37V Z (X) , (45) 

where J is the spin-spin interaction energy and X = tanh(/3J) is the activity of a bond at 
inverse temperature (3. 

(9.) For non-bipartite lattices like the decorated spin lattices there is no direct and fast 
numerical method to compute Pfaffians with their proper sign. This is possible in other 
cases, for example in dimer coverings of bipartite lattices where the matrix M($(a)) is 
block diagonal. In spin lattices we first compute the determinant of M, extract the positive 
root and set Pf ($, 0) positive. The sign at X = 1 is then directly given by cr(a) or obtained 
by analytic continuation of the Pfaffian which is a polynomial of degree 3N in X. Thus cr(a) 
could be used to predict the parity of the number of real zeros of the Pfaffian in the interval 
< X < 1. 



VIII. PRELIMINARY ANALYSIS OF PFAFFIANS 

This section is very preliminary and deals with a number of properties and conjectures 
which are essential in analyzing the behavior of the Pfaffian expansion. Some of these 
properties have been verified on finite Ising lattices and are in agreement with extensive 
numerical sampling and numerical findings ||21|| . We should mention that the procedure has 



not yet been optimized for a numerical approach. In particular the classification of Pfaffian 
symmetries has still to be implemented. 

The self-intersection function cr(a) can be defined on any triangulated surface of genus 
g and particular examples of genus g = to 3 have been worked out in detail in the 
literature |4],|13|]. The argument a takes 2 2g values parametrized by 2g binary variables 
ei, C2 g ■ It is always possible to redefine the basis in Z 2 ) in such a way as to have 

a (a) = (-l)Sti e * e *+s = J]f = i (-l) e ' e >+f. The factor (-l) e »^+9 takes 3 times the value 
1 and once the value —1 as ei,ei +g run on 0,1 and all the factors appearing in a (a) are 
independent. Suppose now that N + (g), N_(g) are the number of times a (a) takes the value 
1,-1 respectively so that iV + (l) = 3 and AL(1) = 1. We have the recursion relation 
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N+(g + l) = 3N + (g) + N-{g), N^g + 1) = N+(g) + 3iV_( 5 ) , (46) 
which has the solution 

N+(g) = 2°- 1 (2° + l), AL(#) = 2^(2* - 1) . (47) 



N + (g) and N_(g) give the total number of positive and negative c(a) in the expansion 
This statistics is important in evaluating the convergence properties of the expansion. For 
high g this means that if we pair off positive and negative cr(a) we are left with a small 
excess of 2 9 positive values, i.e. one part in 2 9 . Since Pf($,0) = 1 we have Zo(0) = 1 as 
expected at T = oo. As X increases all terms in (|44|) become eventually all positive and 
equal to Pf($i?, 1). Comparing the total sum 2 29 Pf($i?, 1) with the known T — > limit for 
Zo we obtain Pf ($f? 1) = 2 14M . The equality of the absolute values of Pf ($, 1) hold only on 
decorated spin lattices whereas generic dimer lattices have a spectrum of values. The factor 
cr(a) induces a cancellation with a cutoff factor 2~ 9 in ( f44"|) in the T — > oo or X — » limit 
which does not occur at T — > and leads to a steeper log derivative for Z (X) as compared 
to that of the single Pfaffians. 

As first noted by Wannier & Kramer |TP|] planar spin lattices can be characterized by 



duality relations. However for g > duality does not relate directly the partition function of 
a lattice to that of the dual but rather acts very simply on the single terms in the expansion 
( [14] ) by changing their signs only so that they are still positive in the T —>■ oo limit. Since 
duality swaps T — ► oo and T — ► limits we see that in dual lattices positivity is required at 
opposite ends of the interval < X < 1. In the 2D Ising lattice the sign reversal does not 
alter the partition function in the thermodynamical limit but the same need not to be true 
in 3 dimensions. In any case the 3D lattice is not self dual. 

A lower bound X for the zeros of Pf ($, X) can be derived from (0). From the spectrum 
of A f ,B f we get \\A f \\ < 2 + V3, \\B f \\ < 2 + V% whence X > ^-^ = 2-^3 ~ 0.267.... 
Numerical analysis indicates that in fact Pfaffians vanish in the range X m = 0.3178... which is 
actually reached by Pf ($^, X) to X M = 0.3506... reached by Pf ($^4, X). These values can 
be computed exactly as roots of an algebraic equation because the translational symmetry 
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of these Pfaffians leads to explicit formulas . Random sampling of orientations $ up to 
Ni — N2 — N 3 — 8 indicate that zeros tend to accumulate around 0.34, see Fig. 5. 
Therefore, we divide now the interval < X < 1 into 4 regions : 



(I) The very low temperature region X > Xm where all terms in the expansion (f4*4|) of 



Zo are positive and the series converges rapidly and agrees with low T expansions. 

(II) The crossover region centered in X = 0.34 where N_(g) = 2 g - 1 (2 g -l) Pfaffians change 
sign. 



(Ill) The region between the estimated critical temperature X = X ~ 0.2108 and X m . 



In this region the negative terms have absolute values smaller that the positive terms 
and we expect that convergence degrades rapidly as we move toward Xq. 

(IV) High temperature region X < X where all absolute values of the terms in ([0]) becomes 
comparable and the value of Z is determined by 2 9 unpaired Pfaffians. For small X 
all absolute values of the Pfaffians are close to 1 and the sum (pj4|) over a sample of L 
random terms has a noise of the order of 2~ 9 a/L while we expect a signal of the order of 
2~ 2g L, normalized to 1 at X = for the complete sum. In order to get a signal we need 
a ratio signal/ noise 2 29 £ ~ 1 i.e. L ~ 2 2g . This means that unless one sums over all 
terms we get only noise and that numerical computation is ruled out unless one obtains 
an explicitly summable formula for the Pfaffians as it happens in the limit T — > 00. 
Therefore exact matching with known T — > 00 results is still possible and useful. As X 
increases the absolute values of positive terms grows on the average more rapidly than 
that of the negative ones thus improving the signal/noise ratio. We conjecture that the 
critical X = X$ is effectively a threshold beyond which the signal becomes effective. 
In general for a fixed X and all $ we have Pf($ A ,X) < Pf($,X) < Pf($ F ,X). 

Details and high/low temperature expansion will be discussed in a forthcoming paper 
2§. 
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IX. DIMER STATISTICS 



A simple application of the above formalism is the evaluations of the number of Perfect 
Matching, i.e. dimer coverings, over the 3D cubic lattice A. 

As for the Ising model, it can be solved exactly or can be treated easily in the case of 
planar lattices whereas it still represents an open problem in the case of non-planar graphs 



The dimer covering generating function is given by a Pfaffian expansion similar to 
where now M($) is a A x A matrix of elements M 9ig /($) = <fii{q) = ±1 depending on the 
orientation of the bond. The decorating bonds are absent whereas the orientations $ of S 
play exactly the same role as in the Ising case. Following the same steps discussed for the 
Ising case and separating odd and even sites, we arrive to a block diagonal form of M($) 



M($) 



( C($) ^ 

-C($) T 



(48) 



Now we have directly Pf(M($)) = Det(C(§)), and the expansion reads: 

ZDimers = ^ <K a ) Pf ( $ ( a )) 

4 E (-l)ELE^^^]^ Det(C(f)), $ = e fc u; fc ) (49) 

{e fc =0,l} k=l 

Such a formula can be used both for exhaustive enumerations of coverings h(N\, N 2 , A 3 ) 
of finite lattices of linear size N% x N 2 X A3 as well as in a probabilistic framework [^6| . 

We have applied (|49[) to the case of finite cubic lattices with open boundaries in order to 
recover and improve the known results. The limitations in the size arise from the number of 
terms appearing in the expansion which increase exponentially with the genus of the surface, 
which for open boundaries grows as g = M 1 M 2 (M 3 - 1) + M 2 M 3 (M 1 - 1) + M^M l (M 2 - 
1) — M\M 2 M% + 1 (g = 2L 3 — 3L 2 + 1 in the isotropic case of linear size L). We have 
found h(4,4,4) = 5051532105 (in agreement with |24j), /i(6,4,4) = 932814464901633 and 
/i(6,6,4) = 123115692449982216049513. 
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Note that the rigorous lower bound [23] to the number of dimer coverings in 3D in the 
M — > oo limit can be easily recovered in our approach by computing, via Fourier Transform, 
the periodic Pfaffian corresponding to $^ [p6f] . 

Eq. (|49| ) can be also thought of as the expansion in terms of determinants of the per- 
manent of — 1 matrices, a j^P — complete problem which can be easily mapped onto the 
evaluation of dimer coverings over an associated bipartite lattice f25[ . 



X. CONCLUSION 

In this paper we propose a combinatorial/topological formalism for the study of the Ising 
problem over lattices of arbitrarily high topological genus which generalizes the well-known 
approach of Kasteleyn. The partition function is written as a sum over Pfaffians with a 
topological signature. 

We apply the method to the 3D cubic Ising problem where we have reached a very 
preliminary assessment on the expansion in the high and low temperature ranges. The same 
formalism applies to the perfect matching problem and provides a determinant expansion 
for the permanent of 0-1 matrices. 

Work is in progress on the physical and algorithmic relevance of the method. 
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FIGURES 




FIG. 1. Orientable surface £ of genus g = N/4 + 1 containing all sites and bonds of the 3D 
cubic lattice A. Some examples of the even (E), odd (O) and exceptional (R) cycles are also shown. 
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FIG. 2. Oriented decorating Clusters. Each site (A) is replaced by either Fisher's planar 
decoration (B) or the equivalent complete graph (C). 
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FIG. 4. 3D simple cubic lattice and the Axial Gauge Tree T 
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FIG. 5. Overlap of the complex zeros of 50 Pfaffians corresponding to different randomly chosen 
orientations (N\ = N2 = N3 = 8). Xm and X m (computed in the N — > 00 limit) are the upper 
ond lower bounds for the zeros given by the singularities of the fundamental and antifundamental 
Pfaffians respectively. Xq corresponds to the value of the critical temperature estimated by different 
analytical and numerical methods. 
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